Generalized rotational Hamiltonians from non-linear angular momentum algebras 
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Higgs algebras are used to construct rotational Hamiltonians. The correspondence between the 
spectrum of a triaxial rotor and the spectrum of a cubic Higgs algebra is demonstrated. It is shown 
that a suitable choice of the parameters of the polynomial algebra allows for a precise identification 
of rotational properties. The harmonic limit is obtained by a contraction of the algebra, leading to 
a linear symmetry. 
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I. INTRODUCTION 



Mathematical aspects of the Higgs algebras have been 
studied intensively in the past [J, UE1 . Generally speak- 
ing, in the framework of non-linear algebras, the com- 
mutation relations between the generators are expressed 
as a linear combination of power terms of them, this is 
the case of the sl^ 3 \2,M.) algebra, where the commuta- 
tor of the ladder operators contains up to cubic powers of 
J3. In spite of their relative involved structures, the fact 
that the algebras contain as limiting cases the more famil- 
iar structures of the groups associated to rotations and 
vibrations turns out to be a particularly appealing fea- 
ture for applications to the group classification of nuclear 
Hamiltonians 1,0. Therefore, they may be useful tools 
to classify observables. The search for physical inspired 
Hamiltonians, which may display definite features about 
polynomial algebras is still open. Concrete applications 
of the formalism have been explored more recently in the 
context of schematic models [6| , to complement previous 
mathematical efforts, like the studies of 0,11,0. Among 
the already studied applications of the concept of non- 
linear algebras, the results of Q can be taken as definite 
motivations for our present effort. 

In this work we continue with the study of applications 
of Higgs algebras, by constructing a Hamiltonian which 
reproduces the behavior of an asymmetric rotor and its 
vibrational limit 0]. We shall show that, while the fi- 
nite dimensional representations of non-linear algebras 
are suitable for the description of rotational-like struc- 
tures, the infinite dimensional representations of a Higgs 
algebra contain vibrational- like structures. 

The paper is organized as follows. In Section II we 
present the essentials of the formalism, and introduce 
both the infinite- and finite-dimensional representations 
of sZ( 3 )(2,R). In Section III. A, we present the basic con- 
cepts which we have developed in the construction of a 
triaxial rotor Hamiltonian. In Section III.B we study the 
vibrational limit of this Hamiltonian by contracting a 
particular, s^ 3 ^(2,R), Higgs algebra. Numerical results 
are presented in Section IV. Conclusions are drawn in 
Section V. 



II. POLYNOMIAL ANGULAR MOMENTUM 
ALGEBRAS 



To start with we shall consider polynomial angular mo- 
mentum algebras up to the third power in J3 [l|, |3j . The 
commutation relations of the s^ 3 -'(2,R) algebra are de- 
fined by [H,i,i 



[J 3 ,J±] = ±J±, 

[J+, J_] = a Jl + P4+ 1 .h + 5, 



(1) 



where the parameters a,0,j, and 8 take arbitrary real 
values. The corresponding Casimir operator reads 



c = J+J- + -4 
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Notice that sl^ (2,R) comprises, as particular cases, 
several well-known algebras. In the linear case (a = = 
0) we recover the following Lie algebras: i) s/(2,R) if 
5 = and 7 takes any non-zero real value; (ii) gl(2) if 
both 7 and 8 take any non-zero real value; (iii) the oscil- 
lator algebra if 7 = and 5 takes any non-zero real 
value; and (iv) the Poincare algebra for 7 = 8 = Q In 
all these cases, when a parameter (either 7 or 6) is differ- 
ent from zero, it can be reduced to ±1 by re-scaling the 
Lie generators. On the other hand, the quadratic alge- 
bra sl^(2, R) is obtained when a = and ^ 0. In this 
case, the linear term 7 J3 can be eliminated by making use 
of the automorphism J3 — * J3 — ip, where ip is a constant 
that has to be fitted. We recall that sl^(2, R) is, for ap- 
propriate values of the constants 0, 7 and 8, the dynam- 
ical algebra for the Dicke and Da Providencia-Schiitte 
models (see and references therein) . We also remark 
that in the cubic case given by a 7^ 0, the same kind 
of automorphism allows us to eliminate the quadratic 
term /3Jf. Moreover, the case with = 8 = corre- 
sponds to the well-known Higgs algebra L l]. This cubic 
algebra was found in the study of the Kepler potential 
on a two-dimensional space with constant curvature and 



2 



has motivated further developments of the representation 
theory of polynomial algebras 0, Q . Finally, note that 
from the complementary viewpoint, the structure con- 
stants a,/3,7,<5 in ([1]) can be interpreted as deformation 
parameters that generate a wide class of structures by 
taking the Poincare algebra (with a — = j = 8 = 0) 
as the origin of the space of deformations. 

As far as the finite-dimensional representations of 
s^ 3 )(2,M) is concerned, these are defined on the (2j + l)- 
dimensional basis spanned by the usual vectors \jm) with 
j = 0, i, 1, §, . . . and m = -j, . . . ,0, . . . , j, like, e.g. Q 

Mi m ) = I — + c j \jm), 

J+\3 m ) = f{m)\jm + q), 

J-\jm) = g(m)\jm - q), (3) 

where q is a positive integer number and c is a real num- 
ber. For a given j the product of functions f(m) and 
g(m) fulfils 

f(m - q)g(m) - f{m)g(m + q) 

(4) 

for m = —j, . . . , j, and the functions / and g also obey 
the constraints 

/(i) = /(i-i) = --- = /(i-? + i) = o, ,,s 
g(-j) = g{-i + 1) = . . . = g{-j + g - 1) = o. [0) 

Some other finite-dimensional representations can be 
found in 

Infinite-dimensional representations of sl t - 3 \2,R) can 
be built on the Hilbert space of number states spanned 
by the set namely: 

J 3 \n) = + \n), 

J+\n) — T(n)\n + q), 

J-\n) =G(n)\n-q). (6) 

The functions J-{n) and Q(n) of Eq. ((6| satisfy the con- 
ditions 

F{n-q)Q{n)-F{n)g{n + q) 

g(o) = g(i)... = g( q -i) = o, 

which can be compared with the conditions given by 
Eqs.©-©. 

The expressions appearing in suggest that the gen- 
erators J± play the role of ladder operators with step q 
acting on the state \jm). Likewise, the results shown in 
equations ^ relate J± to ladder operators of a Lie alge- 
bra to power q [16] . In view of the correspondence of Eqs. 



dS])-© and Eqs. ©-(HI) we stress the fact that the struc- 
ture of s^ 3 ^(2,R) allows for a generalization of the defi- 
nition of the usual angular momentum operators in the 
rotational limit (finite-dimensional representation) and 
also a generalization of the harmonic oscillator in the vi- 
brational limit (infinite-dimensional representations). In 
the next section we shall advanced on these concepts. 



III. GENERALISED HAMILTONIANS. 

If sl^ 3 ' (2, R) is interpreted as a dynamical algebra, then 
the most general Hamiltonian based on such non-linear 
symmetry takes the form 

4 

n = J2 C l J l 3 + C+J+ + C-J-, (8) 

i=0 

where the coefficients Ci as well as C± are arbitrary real 
numbers. 

In order to obtain explicit models, let us introduce the 
following finite-dimensional representation ©-((11) with 
q = 2. The functions /(m) and g(m) are, for this case 

/(m) = ay/ (j + m + TjQ + rn + 2)(j - m)(j - m - 1), 
g(rn) = ay/ (j + m)(j + m- l)(j - m + - m + 2), 

(9) 

where a is a real constant. Notice that these functions re- 
semble the eigenvalues of the square of the linear sl(2, M) 
ladder operators. Next, according to Eq. ©, the de- 
formation parameters a, /?, 7, and <5 in © turn out to 
be 

a = -64a 2 , (3 = 192a 2 c, 

7 = -8a 2 (l + 24c 2 - 2j(j + 1)) , (10) 

S = 8a 2 c(l + 8c 2 - 2j(j + 1)) . 

In this representation, the eigenvalue of the Casimir op- 
erator © reads 

(jm\e\jm) = a 2 + 1) - 2c - 4c 2 ) 

x (j(j + l)-2-6c-4c 2 ). (11) 

The non-vanishing matrix elements of the Hamiltonian 
© are given by 

4 

(jm\n\jm) = J2Ci{^ +tf 

(jm + 2\n\jm)~=C + f(m), 
(jm — 2\H,\jm) = C_ g(jn). 

The above expressions demonstrate that finite- 
dimensional representations of non-symmetric algebras 
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are, therefore, suitable for the treatment of rotational- 
like structures. This property is a result of the explicit 
dependence of the functions / and g upon j, which is 
acting as the effective level of the representation since 
f(m) and g(m) vanish if m = j, j ' — 1 and m = —j, —j + 1, 
respectively. 

Similarly, for the case of infinite-dimensional represen- 
tations one has 



the intrinsic angular momentum components respect to 
the body-fixed frame, satisfy the commutation relations 
— —isijklk- We define the angular momentum 
rising and lowering ladder operators 

I ± =h±ih, (19) 
which obey the commutation rules 



T{n) = b^/(n + l)(n + 2), 
Q(n) = b^/n(n - 1), 



(13) 



provided that q = 2, & is a real constant, and the defor- 
mation parameters are 

a = /3 = 0, 7 = -8& 2 , 5 = 26 2 (4c- 1). (14) 
Hence the eigenvalue of the Casimir operator © is 



(n\C\n) = 26 2 c(l + 2c). 



(15) 



For this case the non-vanishing matrix elements of the 
Hamiltonian JEI) are 



(n\H\n) 



4=0 



(n + 2|W|n) = C+T(n), 
(n-2|Wn) =C-Q(n), 



(16) 



to be compared with (|12j) . Thus, infinite-dimensional 
representations of non-linear algebras exhibit vibrational- 
like structures, under a proper choice of parameters, 
based on the explicit dependence of T and Q upon n. 



A. Triaxial rotor 

The results of the previous subsection illustrate rather 
explicitly the scope of the formal connection we would 
like to establish between finite (infinite) representations 
of a polynomial algebra and the generators of angular 
momentum (vibrational quanta) and the corresponding 
Hamiltonians. 

Let us start with the case of a triaxial rotor [J] , whose 
Hamiltonian is written 



where 



1 

21,' 



(17) 



(18) 



The operators Ii are the components of the angular mo- 
mentum along the principal axes, and Xi are the com- 
ponents of the tensor of inertia. The generators X 



of 



[h,I±]=Tl±, [/+,/_] = -2J 3 . 
The Casimir invariant is given by 

3 



(20) 



1 



Ii + §(/+/- +/-/+)• (21) 



12 

i=l 

In this way the Hamiltonian ((TTJ) can be written as 
H 



^(A 1 +A 2 )I 2 + ^(A 1 -A 2 )(ll + I 2 _) 



+ 1(A 1 -A 2 )kII 



where 



2A 3 ~A X -A 2 
A x - A 2 



(22) 



(23) 



is the factor which measures the asymmetry of the rotor. 

In the following, we shall show that a correspondence 
between the triaxial rotor case and a cubic Higgs alge- 
bra can indeed be established by a proper choice of the 
parameters of the non-linear symmetry. 

Firstly, we adopt the following choice of parameters, 
in Eq.@, 

a = -^. = 0, 7=^(2/ 2 + 2/-l), <5 = 0, 

(24) 

for the non- linear symmetry, and taking q = - 2, c = 0, 
and a — jj in Eq.(f4]) one obtains 

f(M) = ay/(J + M+ 1)(J + M + 2)(J - M)(J - M - 1) 
g(M) = ay/(J + M){J + M - 1)(J - M + 1)(J - M + 2), 

(25) 

where J — I and M — Mj, being / the total angular 
momentum and Mj its projection. The operators J± 
and J3 are written 



■> — —I±i J3 — — — ia 



41 2 

The corresponding commutation relations are 

[J 3 ,J±] = ±J±, 

[J+, J-] = (2/ 2 + 27 - 1) J 3 - 4 Ji 



(26) 



P 



(27) 
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Secondly, the Hamiltonian of Eq. (|22|) is expressed in 
terms of the non-linear symmetry as 



1 



H = -{A 1 + A 2 )J{J + 1) + J{A 1 -A 2 ){J+ 

+2(A 1 - A 2 )U$. 

The matrix elements of the Hamiltonian of Eq. 
given by 



J-) 
(28) 



1 



1 



(J M|W|JM> = -{A x + A 2 )J{J + 1) + -(Ai - A 2 )kM , 
(J M + 2|W|JM) = J{A 1 - A 2 )/(M), 
(J M - 2|W|JM) = J{A X - A 2 )g{M). 

(29) 

Both Hamiltonians, the one of Eq. ([22]) and the one of 
Eq. (J2HJ) 7 have the same matrix elements, and in conse- 
quence the same spectrum and the same eigenfunctions. 

The above example shows that the standard angular 
momentum operators of the triaxial rotor can be under- 
stood as the mapped version of the generators of a partic- 
ular non-linear symmetry. That is to say that the triaxial 
rotor is a particular case of a finite-dimensional represen- 
tation of sl^(2,R). 



B. Harmonic limit of the triaxial rotor. 

We shall now illustrate the use of the concept of con- 
traction of a Higgs algebra, in connection with the struc- 
ture of H, of Eq. ([2"2")) . in the case of large angular mo- 
mentum I. In the standard representation of angular 
momentum [Tl| . the Hamiltonian of equation (|22p can 
be analyzed by making use of boson mapping techniques 

We shall use the well known Holstein-Primakoff boson 
mapping to get the images of the operators of Eq. (|2"0|) , 



C f vg - C C 

/_ = y/2I-C({, 

h = i -Pe- 



rn 



where (£) are boson creation (annihilation) operators. 
In the large angular momentum limit, I » 1, the Hamil- 
tonian of Eq. (|2"2")1 reads 



H 



— Hqq + Hu + H 2a + H nn + H , 
= A 3 1(1+1) -^(A 1 -A 2 ) K I, 



--(A x -A 2 )k{2I-1) c f C, 



H 20 = Z( Al -A 3 )I( 



H n 



H' 



(A 1 - A 2 )k [CO 



8/ 



(Ai - A 2 ){ C t3 C + C f C 3 )- (31) 



The terms of H of Eq. (fBTJ) are classified in powers of J, 
and by the corresponding powers of the boson creation 
and annihilation operators p and £ 



(32) 



kll 



The leading order terms, -STrpa = #00 + Hn + H 2 q, can 
be diagonalized by adopting the Random Phase Approx- 
imation (RPA) formalism [H, [l2j ■ This is achieved by 
introducing the RPA boson operator 



rt = x c f - y c, 

and by solving the equation of motion 



(33) 



(34) 



Details of the procedure, which is rather well known, can 
found in textbooks [3, [l3[ • In this approximation Hrpa 
takes the form 



H KPA = A 3 I(I + 1) 



r t r . 



with 



21 y/ (A 3 - Ai)(A 3 - A 2 ). 



(35) 



(36) 



The expression for the RPA eigenfrequency oj is obtained 
by expressing -ffRPA and Pf in terms of the operators C 
and £ and by equaling terms after performing the com- 
mutation (|3~4")> . 

The remaining terms H nn and H ' of the Hamiltonian 
of Eq. (f3Tj) . can then be transformed to the RPA bo- 
son basis (r^, r). By keeping terms up to (r^T) and 
neglecting negative powers of the angular momentum in 
the coefficients ak n {I), one gets 

4,1(1 + 1) + L + i-(A!-A 2 )) (rtr + 0, 
+ ^(Ax-A^^rtr) 2 . (37) 

Proceeding now in a manner analogous to the previ- 
ous subsection, we shall show that this Hamiltonian is 
obtained as a particular case of the sl^(2,R) algebra. 
As we shall demonstrate, the harmonic limit of the tri- 
axial rotor can alternatively be obtained through a con- 
traction procedure and next embedded, as a particular 
case, within the s^ 3 - ) (2,R) structure together with some 
infinite-dimensional representations. In what follows we 
describe explicitly this procedure and, furthermore, we 
shall show that this is equivalent to the RPA approxima- 
tion. It means that, at leading order, both constructions 
lead to the same results. 



5 



The contraction of the Lie algebra ([20|) to the harmonic 
oscillator algebra /14 is obtained by introducing the new 
generators defined by [Ti ]: 

f± = el±, U = -h + -^, (38) 

where e is the contraction parameter and 1 is the identity 
operator, which obviously commutes with I± and I3. By 
computing the new Lie brackets and next applying the 
limit e — ► we find the defining relations for ft, 4 : 

[44] = ±4, 4,4] = i= [V] = 0. (39) 

Since the eigenvalue of 1 is always 1 in all the represen- 
tations here considered, hereafter we always fix 1 = 1. In 
this respect, notice that, in fact, such an identity operator 
should also be introduced in the commutation relations 
of sl^ (2, R) fl} multiplying the parameter 8, but due to 
the above reason this is not usually written explicitly in 
the literature [ESS]. 

The contracted Casimir I 2 comes from (l2Tj) as the limit 

t 2 = lim e 2 (V - ^4) = ~h + 1(1 J- + Ll+). (40) 
The contraction of the rotor Hamiltonian l|22p follows as 

= I(A 1 +A a )i a + i(A 1 -A 2 )(i*+£) 

-i(A, - A 2 )kI 3 , (41) 

where k is the same of Eq. (f2"3")l . Now all of these results 
can be introduced in sl^ (2, M) by taking the parameters 
AH) with c = 0, 

a = P = Q, 7 = -8b 2 , 6 = -2b 2 , (42) 

which implies that the eigenvalue of I 2 is 1/2. The gen- 
erators of s^ 3 )(2,R) are then chosen as 

J±=bl 2 ± , h=\h, (43) 

and fulfil the commutation relations 

[J 3 , J±] = ±J±, [J+, J-] = -8b 2 . h - 46 2 I 2 . (44) 

Hence the harmonic limit of the rotor Hamiltonian arises 
as the sZ( 3 )(2,R)-system given by 

U = 1 -{A, + A 2 )I 2 + l(Ai - A 2 )(J + + J_) 

-(A 1 -A 2 )kJ 3 , (45) 

which compares with the Hamiltonian of Eq. (|28|) up to 
the term proportional to J3, due to the contraction (|38|) . 
Except for an additive constant, Eq. (|4"5"]) coincides with 
the terms of Eq. (|3"Tj) labelled as Hrpa- 



IV. RESULTS AND DISCUSSIONS 

In this section we shall present and discuss the results 
of our calculations. We shall begin with the analysis of 
the rotational regime by calculating the eigenvalues of 
the Hamiltonian of Eq. (|28|) and then we shall proceed 
by showing the results corresponding to the eigenvalues 
of the Hamiltonian of Eq. (j45|) . Since we have already 
presented the algebraic details in the previous sections II 
and III, we shall focus on the numerical results. 

The spectrum of a triaxial rotor has remarkable fea- 
tures [ij like the breaking of degeneracies and the re- 
ordering of levels in degenerate pairs. The transition be- 
tween both regimes is governed by the asymmetry pa- 
rameter k. These features are exhibited, too, by the 
spectrum of eigenvalues of the Hamiltonian (J2SJ). Figure 
1 shows the spectrum as a function of the asymmetry n 
and for a fixed angular momentum, / = 10. Figure 2 
shows the low energy region of the spectrum for 7 = 4. 
As it can be seen from both figures, the just mentioned 
features emerge clearly. It means that the choice ([24]) 
of the non-linear algebra sl^(2,R) leads to the triaxial 
rotor. This is an expected result, because using the al- 
gebraic procedure of Section III we have shown that the 
matrix elements of both Hamiltonians do indeed coincide. 
Figure 3 shows the eigenvalues of the Hamiltonian of Eq. 
(|45[) as a function of the asymmetry k. The spectrum is 
clearly harmonic for a fixed value of k. The dependence 
of the frequency to upon k and angular momentum / is 
a unique feature of the large angular momentum limit of 
the triaxial rotor. The domains \k\ < 1, of the rotor, and 
|k| > 1, of the vibrator are not connected. 

V. CONCLUSIONS 

In this work we have studied the correspondence be- 
tween a polynomial algebra, of the cubic Higgs-type, and 
the standard angular momentum and harmonic oscilla- 
tor algebras. Particularly, we have shown that: (i) by 
a suitable choice of the parameters entering the defini- 
tion of the cubic algebra, one may obtain the triaxial 
rotor, and, (ii) in the limit of large angular momentum, 
the vibrational structure emerges as a contraction of the 
algebra. 

We think that these results suggest the existence of 
direct physical realizations of the Higgs algebra for the 
case of Hamiltonians describing nuclear vibrations and 
rotations, as it was found for the case of classical Hamil- 
tonians, e.g. the study of the Kepler potential on a two- 
dimensional curved space. 

For the specific application to nuclear structure prob- 
lems, the potential of non-linear algebras in dealing with 
the construction of Hamiltonians becomes evident. The 
use of non-linear algebras, like the s^ 3 ^(2,R) or even al- 
gebras with a larger number of generators, may shead 
light on more involved structures lying in between the 
rotational and vibrational extremes. This may be the 
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case of the newly discovered symmetry reported in [1 51 ] - 
Work is in progress concerning this point. 
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FIG. 1: Eigenvalues of the Hamiltonian of Eq. (J28J) , as a func- 
tion of the asymmetry n. The values are shown in arbitrary 
units. The parameters used in the calculations are: / = 10, 
A\ = 1.66, and A2 = 0.83. The results coincide with the ones 
of a triaxial rotor [J 



[16] Generally speaking a boson-like of the form J+ cc a 
and J- oc a q . 
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FIG. 2: Idem as Figure 1, for 1 = 4. 
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FIG. 3: Eigenvalues of the Hamiltonian of Eq, (|45p . as a func- 
tion of the asymmetry k. The parameters are the same given 
in the captions to Figure 1. In curves are shown the results 
of the harmonic limit of Eq. (I28|l obtained by contracting the 
sZ (3) (2,R) algebra. 



